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ABSTRACT 

Phase-resolved VLT spectroscopy of the bursting Low Mass X-ray Binaries 4U 1636- 
536/V801 Ara and 4U 1735-444/V926 Sco is presented. Doppler images of the Mil A4640 
Bowen transition reveal compact spots which we attribute to fluorescent emission from the 
donor star and enable us to define a new set of spectroscopic ephemerides. We measure 
K cm = 277 ± 22 km s" 1 and K cm = 226 ± 22 km s^ 1 from the Mil spots in V801 Ara 
and V926 Sco respectively which represent strict lower limits to the radial velocity semi- 
amplitude of the donor stars. Our new ephemerides provide confirmation that lightcurve max- 
ima in V801 Ara and likely V926 Sco occur at superior conjunction of the donor star and 
hence photometric modulation is caused by the visibility of the X-ray heated donor. The ve- 
locities of Hell A4686 and the broad Bowen blend are strongly modulated with the orbital 
period, with phasing supporting emission dominated by the disc bulge. In addition, a reanaly- 
sis of burst oscillations in V801 Ara, using our spectroscopic To, leads to K\ = 90 — 113 km 
s _1 . We also estimate the K — corrections for all possible disc flaring angles and present the 
first dynamical constraints on the masses of these X-ray bursters. These are K2 = 360 ± 74 
km a -1 , / (M) = 0.76 ± 0.47 M and q = 0.21 - 0.34 for V801 Ara and K 2 = 298 ± 83 
km s" 1 , f (M) = 0.53 ± 0.44 M Q and q = 0.05 - 0.41 for V926 Sco. Disc flaring angles 
a > 12° and q ~ 0.26 — 0.34 are favoured for V801 Ara whereas the lack of K\ constraint 
for V926 Sco prevents tight constraints on this system. Although both binaries seem to have 
intermediate inclinations, the larger equivalent width of the narrrow NIII line in V801 Ara at 
phase 0.5 relative to phase suggests that it has the higher inclination of the two. 

Key words: stars: accretion, accretion discs - binaries: close - stars: individual: (V801 Ara; 
V926 Sco)- X-rays: binaries - 



1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting binaries where 
a low mass donor transfers matter onto a neutron star or black 
hole. 4U 1636-536 (=V801 Ara) and 4U 1735-444 (=V926 Sco) 
are among the optically brighter members in the class of persistent 
LMXBs, characterized by L x ~ io 37-38 erg s _1 and blue spectra 
with weak high-excitation emission lines (mainly Hell A4686 and 
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the CIII/NIII Bowen blend at A4640). They also share similar prop- 
erties: they are both atoll sources ( as based o n the patt ern traced in 
X-ray color-color diagrams, see e.g lHasinger & van derKlisI 19891) 
with frequent burst activity and short orbital periods (3.80 and 4.65 
hrs respectively) revealed th rough optical photometry ICorbet et alJ 
1986, Pede rsen et alj |1981). Their lightcurves display shallow si- 
nusoidal modulations which have been interpreted as due to the 
geometrically varying visibility of the irradiated donor star (e.g. 
Ivan Paradiis et alJll988t) . Therefore, the photometric maxima sup- 
posedly define orbital phase 0.5 i.e. inferior conjunction of the 
compact object. Note, however, that this assumption requires con- 
firmation because photometric maxima can sometimes be associ- 
ated with asymmetries in the disc structure such as the visibility of 



2 J. Casares et al. 



the irradiated inner disc bulge at phase ~ 0.3 (e.g. 4U 1822-371 
|H^eTjig r & Masonl|l989l) or superhump activity (see lHaswell et al] 

Only a few spectroscopic studies have been presente d on th ese 
two binaries up to now. For example, ISmale & Corbetl 1 1 99 ll) re- 
port He spectroscopy of V926 Sco showing that the line core 
is dominated by emission from the disc bulge or splash region 
where the gas stream interacts with the outer disc rim. On the other 
hand. lAugusteiin et alJil99^) (A98 hereafter) present radial veloc- 
ity curves of Hell A4686 and the Bowen blend at A4640 for both 
V926 Sco and V801 Ara. They conclude that these high excitation 
lines are also tracing the motion of the disc bulge. 

V801 Ara is particularly remarkable since it is one of only 
14 bursters where "burst oscillations" (i.e. nearly coherent high- 
frequency pulsatio ns) have been dete cted during several thermonu- 
clear X-ray bursts iGiles et alJ2"ool G02 hereafter). Furthermore, 
a train of burst oscillations was also discovered in a 13 min in- 
terval during a "superburst", showing a frequency drift which 
has been interpreted as the orb ital motion of the neutron star 
I Stroh maver & Markwardtl f2Q02. SM02 hereafter). By fitting the 
frequency evolution with a circular orbit SM02 constrain the radial 
velocity amplitude of the neutron star to the range 90 < K\ < 175 
km s _1 . These constraints on K\ will be readdressed and improved 
in this paper. Further constraints on the system parameters and the 
component masses, however, require dynamical information on the 
donor star which, unfortunately, is normally overwhelmed by the 
optical emission from the X-ray irradiated disc. 

A new indirect method to extract dynamical informa- 
tion from donor stars in persistent LMXBs was proposed by 
ISteeehs & CasareH J20 02) (SC02 hereafter). They detected narrow 
high excitation emission lines from the irradiated surface of the 
donor star in Sco X-l which led to the first determination of its 
radial velocity curve and mass function. These lines are strongest 
in the Bowen blend, a combination of C III and N III lines which 
are powered by photoionization and fluorescence (respectively) due 
to UV photons arising from the hot inner disk l McC lintock et alJ 
Il975l) . This technique has been amply confirmed by subsequent 
studies of the eclipsing ADC (Accretion Disc Corona) and X-ray 
pulsar 4U 1822-371 JCasares et all2003l) . t he black hole soft X -ray 
transient GX 339-4 during its 2002 outburst 1 Hvnes et al. 2003) and 
Aql X-l during its 2004 outburst iCorneli sse et alj2006fc 

In this paper we apply this method to the X-ray bursters 4U 
1636-536 (=V801 Ara) and 4U 1735-444 (=V926 Sco) and present 
the first detection of the donor stars in these two LMXBs. The pa- 
per is organized as follows: Section 2 summarizes the observation 
details and data reduction. The average spectra and main emission 
line parameters are presented in Sect. 3, with multigaussian decon- 
volution of the Bowen blend. In Section 4 we analyse the radial 
velocities and orbital variability of the strong Bowen blend and 
Hell A4686 emission lines. Estimates of the systemic velocities 
are obtained through the Double Gaussian technique applied to the 
wings of the Hell line. Using these systemic velocities we compute 
Doppler tomograms of Hell A4686 and the Bowen fluorescence 
NIII A4640 which are presented in Section 5. The NIII maps dis- 
play evidence of irradiated emission from the donor star, which is 
used to refine the absolute phasing and systemic velocities. Finally, 
in Section 6 we provide an improved determination of K\ for V801 
Ara and present our constraints on the masses in the two binaries. 
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Figure 1. Summed spectra of V801 Ara and V926 Sco with the principal 
emission lines indicated. 

2 OBSERVATIONS AND DATA REDUCTION 

V801 Ara and V926 Sco were observed on the nights of 23 and 
25 June 2003 using the FORS2 Spectrograph attached to the 8.2m 
Yepun Telescope (UT4) at the Observatorio Monte Paranal (ESO). 
A total of 42 spectra of 600s of V801 Ara and 102 exposures of 
200s of V926 Sco were obtained with the R1400V holographic 
grating, covering a complete orbital cycle per night for each tar- 
get. We used a 0.7 arcsec slit width which rendered a wavelength 
coverage of AA4514-5815 at 70 km s -1 (FWHM) resolution, as 
measured from Gaussian fits to the arc lines. The seeing was vari- 
able between 0.6" - 1.2" during our run. The flux standard Feige 
110 was also observed with the same instrumental configuration to 
correct for the instrumental response of the detector. 

The images were de-biased and flat-fielded, and the spectra 
subsequently extracted using conventional optimal extraction tech- 
niques in order to optimize the signal-to-noise ratio of the out- 
put spectra (Home 1986). A He+Ne+Hg+Cd comparison lamp im- 
age was obtained in daytime to provide the wavelength calibration 
scale. This was obtained by a 4th-order polynomial fit to 19 lines, 
resulting in a dispersion of 0.64 A pix^ 1 and rms scatter < 0.05 
A. Instrumental flexure was monitored through cross-correlation 
between the sky spectra and was found to be very small, always 
within 14 km s _1 (0.4 pix.) on each night. These velocity drifts 
were removed from each individual spectrum, and the zero point 
of the final wavelength scale was established from the position of 
the strong OI A5577.338 sky line. All the spectra were calibrated 
in flux using observations of the flux standard Feige 110. However, 
due to light loss caused by our narrow slit and variable seeing con- 
ditions, our flux calibration is only accurate to ~ 50%. 



3 AVERAGE SPECTRA AND EMISSION LINE 
PARAMETERS 

Figure Q presents the average spectra of V801 Ara and V926 Sco 
in fx flux units. They show a blue continuum with broad high exci- 
tation emission lines of Hell A4686, A5411, the Bowen blend at 
AA4630-50 and H/3, which are typical of X-ray active LMXBs. 
Possible Hel lines at A4922 and A5015 are also identified but 
these are significantly weaker. Table 1 summarizes the FWHM, EW 
and centroid A c of the main emission lines obtained through sim- 
ple Gaussian fits. We note that the average spectra (including the 
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Table 1. Emission line Parameters 



Line FWHM EW Centroid AV 

(kins' 1 ) (A) (A) (kins' 1 ) 

V801 Ara 

Bowen 1848 ± 65 4.25 ± 0.03 4643.2 ± 0.4 

HeIIA4686 1216 ±45 3.31 ± 0.03 4685.6 ± 0.3 -9 ±11 

H/3 963 ± 56 1.67 ± 0.02 4860.9 ± 0.3 -25 ±21 



V926 Sco 

Bowen 1662 ± 65 3.99 ± 0.02 4641.5 ±0.4 

HeIIA4686 657 ± 26 2.18 ±0.02 4683.9 ± 0.2 -118 ±12 

HP 507 ± 39 0.43 ± 0.02 4859.2 ± 0.4 -134 ± 25 



intensities. The line widths are mainly driven by the fit to the un- 
blended Hell line and hence the latter is effectively used as a tem- 
plate to constrain the line profiles within the Bowen blend. Figure|2| 
presents the results of the fit. Because widths are set to be the same 
for all lines, flux ratios are given by simple peak ratios. Our best 
fit yields NIII ratio J(A4634)//(AA4641 - 2) = 0.31 ± 0.10 and 
0.59 ± 0.02 for V801 Ara and V926 Sco, respectively. This is a 
fa ctor < 2 lower than the theoretical NIII ratio of 0.71, computed 
by Nussbaumer 1 1 97 1 ) for Bowen fluorescence, but consistent with 
the results of Schachte r et aljfT989l) for Sco X-l. Similarly, Hynes 
(private communication) finds a NIII ratio in the range 0.36-0.56 
(depending on time) for GX 339-4 during its 2002 outburst. We 
can also calculate the CIII/NIII ratio /(AA4647 - 52)/I(A4634 + 
AA4641 - 2) and find 0.38 ± 0.06 and 0.35 ± 0.01 for V801 Ara 
and V926 Sco respectively. For comparison, we have performed 
the same analysis on average spectr a of Sco X-l a nd 4U 1822-371, 
using data presented in SC02 and Casares et al. 1 2003), and find 
0.44 ± 0.03 and 0.49 ± 0.01 respectively. The lower CIII/NIII ra- 
tio in V801 Ara and V926 Sco may indicate possible evidence of 
CNO processed material but we have to be cautious here because of 
the oversimplification in our fitting model. In particular, we cannot 
rule out other possible transitions (most likely Oil lines at A4641.8 
and A4649.1; see iMcClintock et all 19751) contributing to our com- 
ponents at ~464lA and ~465lA, which so far we have assumed 
to be dominated by NIII and CIII transitions. Furthermore, since 
the CIII and Oil lines are not powered by Bowen fluorescence but 
photoionization, different line ratios may stem from differences in 
efficiency of these two excitation mechanisms rather than true CNO 
abundance variations. 



4 RADIAL VELOCITY STUDY AND ORBITAL 
VARIABILITY 

Figure [3] displays the radial velocities of the Bowen blend and 
Hell A4686 line for V801 Ara and V926 Sco, obtained by cross- 
correlating the individual spectra with Gaussians of fixed FWHM 
as given in Table 1. We have assumed A c = 4643.0 for the central 
wavelength of the Bowen blend. The velocities are folded in orbital 
phase, using the ephemerides of G02 and A98 for V801 Ara and 
V926 Sco but shifted in phase by +0.5 to make phase zero coin- 
cide with the inferior conjunction of the donor star. These will be 
called the photometric ephemerides and will be used throughout 
this Section. The accumulated phase uncertainty at the time of our 
observations is ±0.06 for both systems. Note that this phase con- 
vention assumes that lightcurve maxima are driven by the visibility 
of the irradiated donor star. 

The radial velocity curves for V801 Ara show maxima at 
phase ~ 1. This phasing suggests that the velocity variations 
contain a significant component arising in the disc bulge (which 
has its maximum visibility at phase ~ 0.75), caused by the in- 
teraction of the gas stream with the outer edge of the disc, 
as is typically seen in persistent L MXBs (e.g. 4U 1822-371; 
ICowlev. Crampton & Hutchingslll982l) . On the other hand, the ra- 
dial velocity curves for V926 Sco show maxima at phase ~ 0.7. By 
comparing the Hell curves in the two binaries we clearly see evi- 
dence for a large system ic velocity of ~ -150 km s _1 in V926 Sco, 
in good agreement with Smale & Corbej 119911) . We also note that 
the CIII/NIII blend in V801 Ara is modulated with a much larger 
amplitude than in V926 Sco i.e. 280 km s _1 versus 70 km s _1 , 
respectively. 

Figure [4] presents the trailed spectra of the Bowen blend and 
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Figure 2. Combined fit to Hell A4686 and the Bowen blend using three 
Gaussians at ~4634A (NIII), ~464lA (NIII) and ~465lA (CIII). See 
text for details. 

line strengths) look very similar to those presented by A98, which 
seems to imply that there has been no large, long-term variations 
between the two data epochs. Incidentally, the EWs and FWHMs 
of all lines do not show any significant nigh-to-night variability 
nor modulation with the orbital period. Aside from the intrinsically 
broad Bowen blend, which is a blend of at least three CIII/NIII 
transitions, we note that emission lines are a factor ~2 narrower 
in V926 Sco than in V801 Ara. Given the similarities in their or- 
bital periods, this suggests a projection effect with a lower inclina- 
tion angle for V926 Sco. In addition, the emission line centroids 
in V926 Sco are significantly blueshifted, probably due to a larger 
(approaching) systemic velocity. This also applies to the Bowen 
complex since the difference in wavelength between V801 Ara and 
V926 Sco (-1 10 km/s) is consistent with the difference in velocities 
found for the other lines. The blue continuum is also steeper for 
V926 Sco, as expected because of its lower reddening (see A98). 

The Bowen blend mainly consists of emissions correspond- 
ing to the NIII tr ansitions at A A4634,4 64 1,4642 and CIII at 
A A4647 ,465 1,4652 fcchachter etailll989t) . In an attempt to esti- 
mate the relative contribution of the different transitions we have 
performed a combined multi-Gaussian fit to the average emission 
profiles of Hell A4686 and the Bowen blend for the two LMXBs. 
The fit consists of four Gaussians, one for the Hell line, two for 
the NIII emissions at A4634 and A4641-2 and another for the CIII 
emissions at A4647-52. Free parameters are the line widths, which 
are set to be equal for all the lines, individual line centroids and 
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Figure 3. Radial velocity curves of the Bowen blend (top) and Hell A4686 
in V801 Ara and V926 Sco. 
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Figure 4. Trailed spectra showing the orbital evolution of the Bowen blend 
and Hell A4686 in 15 phase bins using the photometric ephemerides 

HeIIA4686 for V801 Ara and V926 Sco, after co-adding our indi- 
vidual spectra in 15 phase bins to improve statistics. The core of the 
emission lines in V801 Ara show clear S-wave components. The S- 
wave in the Bowen blend is quite narrow and it may arise from the 
heated face of the donor star, as observed in Sco X-l (SC02). The 
S-wave in Hell is rather extended and is likely produced in a region 
with a large velocity dispersion such as the disc bulge. On the other 
hand, the Bowen blend in V926 Sco is rather noisy and does not ex- 
hibit any clear components that are visible by eye. The Hell trailed 
spectra do show a clear orbital modulation with complex structure. 

Following our work on Sco X-l (SC02) and 4U 1822-371 
ICasares et alJl2003l). we have applied the double-Gaussian tech- 
nique Isrimeider" &. Yound fl98CO to the wings of the Hell A4686 
line in an attempt to estimate the radial velocity curve of the com- 
pact object. The trailed spectra presented in Fig. [4] demonstrate 
that the line cores are dominated by strong, complex, low-velocity 
components associated with asymmetric emission from the outer 
disc and/or the donor star, which we want to avoid. Therefore, 
by convolving the emission line with a double-Gaussian filter of 
sufficiently large Gaussian separation we can extract radial veloc- 
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Figure 5. Diagnostic Diagram for Hell A4686 in V801 Ara (solid circles) 
and V926 Sco (open circles). The dotted horizontal line marks <j>o = 0.5 
or the expected inferior conjunction of the compact object, according to the 
photometric ephemerides. 

ity curves from the wings of the profile, which are expected to 
follow the motion of the compact star. We have used a double- 
Gaussian bandpass with FWHM = 100 km s _1 and relative 
Gaussian separations in the range a = 400- 1400 km s _1 in steps 
of 100 km s _ . In order to improve statistics, we have co-added 
our spectra in 15 phase bins using the photometric ephemerides. 
The radial velocity curves obtained for different Gaussian separa- 
tions are subsequently fitted with sine-waves of the form V(<fi) = 
7 + K sin 27r(0 — <f>o), fixing the period to the orbital value. The 
best fitting parameters are displayed as a function of the Gaussian 
separation a in a Diagnostic Diagram (see Fig.[5}. 

In both cases we see how the line cores are dominated by high- 
amplitude (K ~ 200 km s" 1 ) S-waves which fade as we move to 
the line wings, where lower K— amplitudes of a few tens of km 
s _1 are found. On the other hand, the 7 velocities are very steady 
throughout the profiles, with average values of ~ —30 (V801 Ara) 
and ~ —130 km s _1 (V926 Sco). The blue-to-red crossing phase 
for V926 Sco displays a smooth decreasing trend from the line core 
to the wings, whereas it is rather constant for V801 Ara, with an av- 
erage value of 00 ~ 0.7. We estimate that the velocity points start 
to be corrupted by continuum noise for Gaussian separations larger 
than a ~ 10 00 as indicated by the diagnostic parameter a(K)/K 
(see lShafter et al. 1986). Therefore, we decided to adopt the aver- 
age values for the parameters obtained from the last two separa- 
tions before o(K)/K starts to rise i.e. a = 900 — 1000 km s _1 
for V801 Ara and 1000-1100 km s _1 for V926 Sco. This yields 
Ki = 93 ± 25 km s~\ 7 = -42 ± 4 km s"\ <p = 0.68 ± 0.03 



Irradiated donor in4U 1 636-5 36/V801 Ara and4U 1 735-444/V926 Sco 5 



for V801 Ara and Kx = 96 ± 9 km s _ \ 7 = -121 ± 7 km 
s"\ 0o = 0.20 ± 0.01 for V926 Sco. The errorbars have been 
adjusted to incorporate the scatter between different values within 
our preferred range. Note that 4>o of V801 Ara is delayed by ~ 0.18 
with respect to the inferior conjunction of the compact object (dot- 
ted line in Fig. as predicted by the photometric ephemerides. 
This is classically observed in interacting binaries and interpreted 
as contamination of the line win gs by residual emission from the 
disc-bulge/hot-spot <Marshlll998h . On the other hand, (j> of V926 
Sco leads the inferior conjunction of the compact star by ~ 0.30, 
according to the photometric ephemerides. Such a large shift is un- 
expected and a full discussion is diverted to Sect. 6. Furthermore, 
we also note the very asymmetric distributions obtained in the Hell 
Doppler maps (see next Sect.) which may invalidate the double- 
Gaussian technique. Therefore, the constraints on zero-phase, 7- 
velocity and K\ (for V801 Ara) will be superseded in the following 
sections. 



5 BOWEN FLUORESCENCE FROM THE IRRADIATED 
COMPANION 

Doppler tomography enables us to map the brightness distribution 
of a binary system in velocity space through combining orbital 
phase spectra (see details in Marsh 2001). This is particularly ef- 
fective when dealing with weak emission features which are barely 
detected or embedded by noise in individual spectra, such as our 
narrow Bowen emission lines. In order to compute Doppler im- 
ages of the NIII A4640 fluorescence line in V801 Ara and V926 
Sco we have rectified the individual spectra by subtracting a low 
order spline fit to the continuum regions and rebinned them into 
a uniform velocity scale of 37 km s~ per pixel. Doppler images 
were subsequently computed using the photometric ephemerides 
and 7 = -42 and -121 km s" 1 for V801 Ara and V926 Sco, 
respectively. The maps show compact spots shifted in phase by 
+0.03 and —0.18 with respect to the expected location of the donor 
star in the velocity maps (i.e. along the vertical V y axis) 1 . Assum- 
ing that these spots are produced on the irradiated hemisphere of 
the donor star, as h as been shown to be t he case in Sco X-l (SC02) 
and 4U 1822-371 (Casares et al. 2003J), we correct the previous 
photometric ephemerides and derive the following spectroscopic 
ephemerides, which will be used in the remainder of the paper: 

Tq(HJD) = 2452813.531(2) + 0.15804693(16)S (1) 

Tq(HJD) = 2452813.495(3) + 0.19383351(32)^ (2) 

where the zero-phase error comes from the uncertainty in the 
centroid position of the spots. Equation (1) corresponds to the 
ephemerides of V801 Ara and equation (2) to V926 Sco. With these 
ephemerides, the Hell A4686 maps of the two binaries show a cres- 
cent shape brightness distribution pointing towards emission from 
an extended disc bulge. 

Following SC02 we have used the Hell Doppler maps to re- 
fine the systemic velocities derived in the previous section. The \ 2 
value of the map was calculated for a range of 7's, and the best 
fit in terms of mimimal \ 2 was achieved for 7 ~ — 39 km s -1 
(V801 Ara) and 7 ~ -132 km s^ 1 (V926 Sco). The use of an 
incorrect systemic velocity has the effect of blurring bright spots 

1 Note that the greyscale in the NIII maps have been set to enhance the 
contrast of the sharp components. This is the reason why these maps do not 
show any trace of emission from the underlying broad component. 
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Figure 6. Doppler maps of Hell A4686 (left panels) and NIII A4640 (right 
panels) for V801 Ara and V926 Sco, computed using our spectroscopic 
ephemerides. Systemic velocities 7 = —34 km s" 1 and 7 = —140 km 
s^ 1 were adopted for V801 Ara and V926 Sco, respectively. For compari- 
son we also plot the gas stream trajectory (in units of 0.1 RlI) and Roche 
lobe of the donor star for the case of K2 = 360 km s~ 1 (V801 Ara) and 
K2 = 298 km s~ 1 (V926 Sco). A mass ratio q = 0.27 has been assumed 
for the two binaries. 



into , as opposed to V801 Ara,elongated "defocused" features in 
Doppler images. Therefore, as a further test, we have also com- 
puted NIII A4640 maps for a set of 7-velocities in the range -200 
- +200 km s _1 and we looked for the best focused NIII spots by 
computing the skewness using box sizes of 5, 10 and 20 pixels. We 
find that the most symmetric and compact spots are found for 7 val- 
ues in the range -29 - -39 km s" 1 (V801 Ara) and 137-143 km s _1 
(V926 Sco). Therefore, we decided to adopt 7 = —34 ± 5 km s~ 
(V801Ara) and 7 = -140±3 km s" 1 (V926 Sco) and these are the 
values that we have used in the final Doppler maps, which are pre- 
sented in Fig. [6] The NIII maps reveal compact bright spots which 
we interpret as irradiated emission from the inner hemisphere of 
the donor star. These spots are obviously located along the vertical 
axis because of our ephemerides definition. We have calculated the 
position of these spots by computing their centroids using the cen- 
tering algorithm CENTROID in IRAF 2 and find K em = 277 ± 22 
kms" 1 (V801 Ara)and226±22kms- 1 (V926 Sco). We note that 
K enl is very weakly dependent on 7, with only a 10 km s~ drift 
when 7 varies by ±20 km s" 1 around our central value. There- 
fore, our K em determinations seem very robust and they represent 
strict lower limits to the true radial velocity semiamplitude K2 of 
the donor stars because they must arise from the irradiated hemi- 
sphere facing the neutron star. Further confirmation of our system 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities of Research in As- 
tronomy, Inc., under cooperative agreement with the Nacional Science 
Foundation. 



6 J. Casares et al. 



V801 Ara 



V926 Sco 




niii 




niii | cm 

n I 








1 1 ItttTtttTtttI 





4620 4640 4660 4680 4700 

Wavelength (A) 



4620 4640 4660 4680 4700 

Wavelength (A) 



Figure 7. Average spectra of V801 Ara (left panel) and V926 Sco (right 
panel). Top spectra have been Doppler corrected in the rest frame of the 
donor using the spectroscopic ephemerides and the X-velocity amplitudes 
derived from the bright spots in the NIII A4640 Doppler images. The nar- 
row NIII and CIII lines are clearly seen in the Bowen blend. These are 
smeared out and cannot be detected in the straight averages shown at bot- 
tom. Noise has been filtered using a 2 pixel boxcar smoothing. 




400 600 800 1000 

Time (3) 

Figure 8. Top panel: Dynamic Fourier power contours showing the excur- 
sion in oscillation frequency during the 2001 February 22 superburst from 
4U 1636-53 (V801 Ara) with the best fitting orbit model overlayed. Time 
zero for the x-axi is MJD 51962.7102530038 (TDB) Bottom panel: Pulse 
phase residuals for the best fitting model. 



parameters is provided by the appearance of the sharp CIII/NIII 
transitions after coadding all the spectra in the rest frame of the 
NIII A4640 emission line region (see Fig.0. Note that we have 
not attempted to compute Doppler images of the H/3 lines because 
these are contaminated by phase-variable absorption components, 
as was also the case in 4U 1822-371 (e.g. lCasares et all2003h . Ab- 
sorption violates the principles of Doppler Tomography and makes 
it very difficult to interpret the corresponding Doppler images. 



6 DISCUSSION 

The detection of NIII A4640 fluorescence emission from the irra- 
diated donor in V801 Ara and V926 Sco provides new absolute 
(spectroscopic) ephemerides and opens the possibility to derive the 
first constraints on the dynamical masses of these two LMXBs. In 
the light of the new ephemerides, the photometric lightcurve max- 
ima of V801 Ara take place at phase 0.47 ± 0.06 (with the un- 
certainty dominated by error propagation of the G02 ephemerides) 
and, therefore, they are consistent with X-ray irradiation of the 
donor star. On the other hand, lightcurve maxima in V926 Sco are 
located at orbital phase 0.82 ± 0.06. This is in between the maxi- 
mum visibility of the irradiated donor and the disc-bulge and hence 
it would imply that lightcurve maxima arise by a combination of 
these two emitting regions. However, we note that the uncertainty 
in the photometric ephemerides of V926 Sco is likely to be higher 
than stated in A98 due to the limited number of arrival times fitted 
and their large scatter (Augusteijn, private communication). There- 
fore, it is possible that lightcurve maxima in V926 Sco are also 
consistent with superior conjunction of the irradiated donor star. 
This is the most likely scenario because it would be hard to under- 
stand an emission source peaking at phase ~ 0.8 in an irradiation 
dominated environment without invoking some strange (and stable) 
disc configuration. 

We have also determined the systemic velocities for the two 
binaries i.e. -34 ± 5 km s" 1 (V801 Ara) and -140 ± 3 km s" 1 
(V926 Sco). These can be compared with Vt, the radial velocity, 



relative to the Local Standard of Rest, due to Galactic rotation at 
the location of the binaries. For the case of V801 Ara, we take 
I = 332.9°, b = -4.8° and d = 6 - 7 kpc, based on peak 
fluxes of ra dius-expansion X-ray bursts for Mi = 1.4 — 2 Mq 
(Gallo way et alJl2005l). A ssuming the Galactic rotation curve of 
Naka nishi & SofueH20o3) we find V t in the range - 1 03 to - 1 1 1 km 
s . This is much larger than our systemic velocity which could 
be explained through the recoil velocity gained by the binary dur- 
ing the supernova explosion which formed the neutron star. Thus, 
we can set a lower limit to the kick velocity in V801 Ara of ~ 70 
km s" 1 . Regarding V926 Sco, we take I = 346.1°, b = -7.0° 
and a distance d = 9.1 kpc (for a canonical 1.4 Mq neutron star, 
see A98). The Galactic rotation curve at the position of V926 Sco 
yields V r = —111 km s , close but significantly lower than our 
systemic velocity. The difference might also be ascribed to a kick 
velocity received by the neutron star at birth or, alternatively, a neu- 
tron star mass of ~1.7 Mq. Note that the latter implies a larger 
distance of d = 10 kpc and hence Vt = —140 km s . 

6.1 System Parameters for V801 Ara 

The NIII A4640 spot in the Doppler map yields K em = 277 ± 22 
km s" 1 . It must arise on the inner hemisphere of the donor and 
hence it sets a lower limit to the velocity semi-amplitude of the 
companion's center of mass K-z. In order to find the real K2 
we need to calculate the K — correction or K cul / K2 for the 
case of emission lines in illuminated atmospheres, which depends 
mainly on the binary mass ratio q — M2/M1 (with M2, Mi 
the masses of the donor and compact star respectively) and the 
disc flaring angl e a. The K — correction has b een calculated 
by iMunoz-Daria s. Casares & Martfnez-Paisl 120051) using an irra- 
diation binary code which includes shadowing by an axisymmet- 
ric flared disc. The results are tabulated as a function of q, a and 
the inclination angle i in Table 1 of their paper, although the de- 
pendence on i is very weak. The K — correction is constrained 
between a — (i.e. we neglect disc shadowing) and the geomet- 
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ric limit set by emission from the irradiated li mb of the donor, i.e. 
ijj^/iG < 1 - 0.213<? 2/3 (1 + g) 1/3 (see iMunoz-Darias et alJ 

An upper limit on q is established by taking Ml > 1.4Mq 
and M 2 < 0.48 Mq, the largest possible zero-age mai n sequence 
(ZAMS) star fitting in a 3.79 hr period Roche lobe flout et all 
1996), which leads to q < 0.34. As a first approach, one can as- 
sume the empirical Mass-Radius relation fo r low-mass stars in cat- 
aclysmic va riables and LMXBs (see e.g. ISmifh & Dhillonlll998l 
Warnetl [l995h which yields M 2 — 0.32. This, combined with a 
canonical neutron star of Mi ~ 1.4 Mq, would lead to a plausible 
q ~ 0.23. The K - correction for K em = 277 ± 22 km s" 1 
and q ~ 0.23 would yield 303 < K 2 < 404 km s" 1 , where we 
have a dopted the co efficie nts for the case i = 40° in Table 1 of 
Muno z^Darias et al"H 2005) because V801 Ara is not eclipsing. 

A more refined K — correction requires a knowledge of q. 
This can be constrained by the rotational broadening Kot sin i of 
the companion star which, for the case of synchronous rotation, is 
related to q through 

Kot sini ~ 0.462A- 2 g 1/3 (l + q) 2/3 (3) 

IWade & Hornell988l) . A lower limit to Kot sin i can be estimated 
from the width of the sharp NIII fluorescence emission in the 
Doppler corrected average spectrum. This is because emission lines 
arise not from the entire Roche lobe but from the irradiated part 
only. A multigaussian fit to the Bowen profile presented in Fig.0 
gives FWHM(NIII\4640) = 140±21 km s"\ which includes 
the effect of our intrinsic instrumental resolution (70 km s _1 ). 
This has been accounted for by broadening a Gaussian template of 
FWHM = 70 km s" 1 between 10 and 200 km s" 1 , in steps of 10 
km s" 1 , using a Gray rotational profile iGrav 1992) without limb- 
darkening, because fluorescence lines arise in optically thin condi- 
tions. Our simulation indicates that Kot sin i — 92 ± 16 km s _1 is 
equivalent to FWHM = 140 ± 21 km s" 1 . On the other hand, by 
substituting Kot sini > 76 km s _1 and the K — correction for 
a > 0° into equation (3), one finds a secure lower limit q > 0.08. 

An upper limit to the rotational broadening of the donor star is 
set by Kot sin i — 2irR 2 sin i/P. By assuming that the donor must 
be more evolved than a ZAMS star within a 3.79 hr Roche lobe, 
then R 2 < 0.UR @ <Tout et all 19961) . This, together with i < 78° 
(lack of X-ray eclipses for q > 0.08), yields Kot sini < 138 km 
s" 1 . 

Further constraints are provided by the study of burst oscilla- 
tions which can set limits to the neutron star's projected velocity 
Ki{= qK 2 ). SM02 derived 90 < Ki < 175 km s" 1 by fitting the 
frequency drift of highly coherent X-ray pulsations observed in an 
800s interval during a superburst. They used a circular orbit model 
and fixed the zero phase to the ephemeris of G02. We have used our 
new spectroscopic ephemeris to reanalyse the superburst pulsation 
data in order to better constrain K\. Based on our spectroscopic 
ephemeris, the pulsation interval during the superburst from V801 
Ara comes slightly earlier in orbital phase by 0.032 cycles as com- 
pared to the G02 ephemeris. We fit the pulsation data (see SM02 
for details on the phase timing analysis) to a circular orbit model 
with the reference epoch fixed to our new To, and we also fixed the 
orbital period (using the G02 value). This leaves two free parame- 
ters, the projected neutron star velocity, K\, and the rest-frame spin 
frequency, We find acceptable fits with a reference epoch within 
the ±lcr range for our spectroscopic To. The inferred K\ velocity 
ranges from 90 - 113 km s _1 as To ranges over ±1<j. The best fit 
X 2 values range from 14.6 to 10.6 (with 7 dof) over this same range, 
that is, higher velocities are modestly favored, but given the short- 




Figure 9. Constraints on K 2 and q for V801 Ara. The disc flaring angle 
must be between a = 0° and a = max, whereas the rotational broadening 
is constrained between the observed Kot sin i ~ 76 km s" 1 and 138 km 
s _1 , the maximum allowed for a ZAMS star in a 3.8hr period and for i < 
78° . The shaded area represents the K\ solution obtained from a reanalysis 
of the superburst pulsations. These constrain the binary mass ratio to be 
q = 0.21 — 0.34 and K 2 = 286 — 433 km s _1 . For comparison we 
also mark (dashed line) the solution for a hypothetic disc flaring angle of 
a = 12°. 

ness of the pulse train compared to the orbital period, we do not 
consider these differences as significant. With the reference epoch 
fixed, the statistical error on the velocity is much smaller than the 
range given above, so 90 - 113 km s _1 is a robust range for K\ 
at the ±lcr limits of T . The best fit v ranges from 582.04143 to 
582.09548 Hz, and which may help constrain future searches for 
a persistent millisecond pulsar in V801 Ara. Figure [8]summarizes 
the results of the new timing analysis. 

Fig.[5]summarizes all our restrictions in the K 2 — q parameter 
space i.e. 0° < a < max, Kot sini = 76 — 138 km s _1 and 
Ki = 90 - 113 km s" 1 . The shaded area indicates the region 
allowed by our constraints, which yields q — 0.21 — 0.34 and 
K 2 — 286 — 433 km s _1 . These imply a mass function f(M) = 
Mi sin 3 i/(l + q) 2 = 0.76 ± 0.47 M and hence Mi sin 3 i = 
1.24 ± 0.77 Mq. Clearly the large uncertainties are dominated by 
the error in K em and the uncertainty in the disc flaring angle. 



6.2 System Parameters for V926 Sco 

The same reasoning can be applied to V926 Sco. The spot in the 
NIII A4640 map yields K cm = 226 ± 22 km s" 1 whereas the 
ZAMS Mass-Radius relatio n for a 4.65 hr Roche lobe leads to 
M 2 < O.58M <Tout et alj|l996T) and hence q < 0.41. On the 
other hand, we measure FW H M (N 1 1 1 \4640) = 115 ± 23 km 
s _1 which, after deconvolution of the instrumental resolution using 
Gray rotation profiles as above, yields Kot sin i > 71 ± 21 km s~ . 
This lower limit, combined with the K — correction for a > 0° 
and equation (3), yields q > 0.05 and, hence i < 80°. And, under 
the assumption that the donor star is more evolved than a ZAMS 
star, the Roche lobe geometry implies R 2 < 0.54Rq which, for 
i < 80°, leads to Kot sini < 137 km s" 1 . 

All these restrictions translate into constraints on the K 2 — q 
plane which are presented in Fig. 1101 The allowed region results in 
q = 0.05 - 0.41 and K 2 = 215 - 381 km s"\ depending on the 
value of a. These numbers imply /(M) = 0.53 ± 0.44 Mq and, 
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q 



Figure 10. Same as Fig. |9]but for V926 Sco, but with no pulsation con- 
straints yet. 



Table 2. System Parameters. To indicates zero phase or inferior conjunction 
of the donor star. Orbital periods P or b are from G02 and A98. 



Parameter 


4U 1636-536 


4U 1735-444 




V801 Am 


V926 Sco 


P orb (days) 


0.15804693(16) 


0.19383351(32) 


T (HJD-2452 000) 


813.531 ± 0.002 


813.495 ± 0.003 


7 (km s -1 ) 


-34 ±5 


-140 ±3 


K cnl (km s _1 ) 


277 ± 22 


226 ± 22 


q(M 2 /M 1 ) 


0.21-0.34 


0.05-0.41 


Ki (kms" 1 ) 


90-113 




K 2 (kms- 1 ) 


360 ± 74 


298 ± 83 


f(M) (M ) 


0.76 ± 0.47 


0.53 ± 0.44 


i (deg) 


36-60 


27-60 



thus, Mi sin 3 i = 0.80 ± 0.71 M . Unfortunately we do not have 
any contraints on K\ from burst oscillations and hence our mass 
restrictions are not as well constrained as for V801 Ara. Our best 
estimates of the system parameters for both objects are presented 
in Table 2. 

6.3 Constraints on the inclination 

Further constraints on the stellar masses requires a knowledge of 
the binary inclination. While strict upper limits are set by the ab- 
sence of X-ray eclipses, lower limits can be established by com- 
bining Mi < 3.1 Mo (the maximum mass allowed for a stable 
neutron star, Rhoades & R uffinil 191 4) with our mass function and 
q restrictions. This leads to i — 36° — 74° and i — 27° — 80° 
for V801 Ara and V926 Sc o, respective ly. In addition, the physical 
model of Frank, Kins & Lasota 1 1987) indicates i < 60° due to 
the lack of X-ray dips. 

On the other hand, we mentioned in Section 3 that the fac- 
tor ~2 narrowness of the Hell A4686 profile in V926 Sco with 
respect to V801 Ara indicates a lower inclination. To test this hy- 
pothesis we have produced Doppler corrected averages for V801 
Ara and V926 Sco by coadding all the spectra in two bins cen- 
tered at orbital phases 0.0 and 0.5, when the visibility of the irra- 
diated face of the donor is minimum and maximum, respectively. 
These spectra are presented in Fig. II H and show a marked differ- 



ence for the two binaries: the narrow CIII/NIII lines become signif- 
icantly enhanced around phase 0.5 (top spectra) for V801 Ara, but 
not much difference is seen in V926 Sco between phase 0.5 and 
phase 0. This clearly supports a higher inclination angle in V801 
Ara. Furthermore, the relative contibution of the NIII A4640 line 
with respect to the broad base, as estimated through a multigaus- 
sian fit, is a factor ~ 2 larger for V801 Ara than for V926 Sco. 
This can be taken as the relative contributions of the heated donor 
and disc, and hence, strongly suggests that V801 Ara is seen at a 
higher inclination angle than V926 Sco. This seems at odds with 
the fact that the optical lightcurves in V80 1 Ara and V926 Sco 
display similar amplitudes, A ~ 0.2 m ags Iva n Amero ngen et all 
Il987llvan Paradiis et alll99(j) . However. lde Jong et alJJl996i) have 
shown that lightcurve amplitudes are more sensitive to a than i, 
which might simply imply a thicker disc in V926 Sco and hence 
lightcurve amplitudes cannot be taken as a simple indication of the 
inclination angle. 

6.4 Mass Estimate and Evolutionary Scenarios 

Mever & Mever-Hofmeister <1982h analyzed the effects of X-ray 
heating in the vertical structure of discs and demonstrated that these 
are strongly correlated. Even in the absence of irradiation, they find 
a disc opening angle of ~ 6°. De Jong et al. (1996) modelled the 
effects of irradiation in optical lightcurves and, by comparing with 
observations of LMXBS, they infer a mean disc opening angle of 
12°. Most of their systems are short period LMXBs, just like V801 
Ara and V926 Sco, and hence it seems justified to speculate what 
our system parameters would be in this particular case. 

Regarding V801 Ara, a = 12° leads to q ~ 0.25 - 0.30, 
Ki — 356 — 375 km s _1 (see dashed line in Fig. |9j and hence 
/(M) ~ 0.81 ± 0.07 M and Mi sin 3 i ~ 1.32 ± 0.13 M . 
Asssuming that the donor star is a 0.48 M ZAMS, then M\ ~ 
1.6 — 1.9 M and i ~ 60 — 71°. On the other hand, a canonical 1.4 
M would imply an evolved donor with M2 — 0.35 — 0.42 M 
and i > 71°, which is difficult to reconcile with the absence of X- 
ray eclipses, dips and the low amplitude of the optical modulation 
(see Sect. 6.3). Therefore, a flaring angle of a ~ 12° seems to sup- 
port scenarios with ZAMS donors and massive neutron stars. These 
masses can be reproduced by recent evolutionary scenarios where 
some LMXBs are suposed to descend from binaries with high mass 
donors fechenker & King|2002l) . These are significantly evolved at 
the begining of the mass transfer phase and mass is transferred in a 
thermal timescale until mass ratio reverses, which can result in the 
accretion of several tenths solar masses by the neutron star. This 
wou ld make V801 Ara similar to other LMXBs such as 4U 1822- 
371 iMunoz-Darias et all2005l) . 

However, it is also possible to accommodate a canonical neu- 
tron star with a MS donor and a plausible inclination angle by push- 
ing the disc flaring angle to higher values. For instance, a — 16° 
leads to q ~ 0.30 - 0.34, f(M) ~ 0.51 M and Mi sin 3 i ~ 0.88 
M . Then, for Mi = 1.4 M , M 2 ~ 0.45 M and i ~ 57°. 
These masses would be consistent with standard evolutionary pic- 
tures where donors in LMXBs descend from slightly evolved low- 
mass stars through dynamically stable mass transfer and result in 
slightly undermassive stars. 

Note also that a < 12° are virtually ruled out because this 
would imply higher K2 and f(M) values which result in massive 
neutron stars and very high inclinations. For instance, a = 8° leads 
to f(M) ~ 0.96 M and Mi sin 3 i ~ 1.5 M . Then, for i < 76° 
(lack of X-ray eclipses), Mi > 1.7 M whereas a plaussible i ~ 
57° yields Mi ~ 2.6 M . 
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Figure 11. Doppler corrected averages of V801 Ara (left panel) and V926 
Sco (right panel) in the rest frame of the donor star. Top spectra present av- 
erages between orbital phases 0.25-0.75 and bottom spectra between 0.75- 
0.25. The sharp CIII/NIII components in V801 Ara are clearly more intense 
around phase 0.5 than 0, suggesting a higher inclination binary. Spectra 
have been smoothed by a 2 pixel boxcar. 



The case of V926 Sco is unfortunately unconstrained because 
of the lack of an X-ray mass function. For example, assuming 

a = 12° leads to q ~ 0.09 0.38, f(M) ~ 0.22 0.68 

M and Mi sin 3 i ~ 0.29 - -1.07 M Q . The allowed range in Mi 
and M2 is too wide as to impose any useful restriction on possible 
evoultionary scenarios. More, higher resolution data are required to 
measure the NIII A4640 flux and V I0 t sin i as a function of orbital 
phase, from which tighter constraints on the inclination and disc 
flaring angle can be set. This, together with the determination of the 
X-ray mass function for V926 Sco through pulse delays of burst os- 
cillations and smaller errors in the K arn determinations, is expected 
to provide stronger limits on the stellar masses and the evolutionary 
models for these two LMXBs (see e.g. lMunoz-Darias et alj20o3) . 



7 CONCLUSIONS 

The main results of the paper are summarized as follows: 

• We have presented the first detection of the donor stars in the 
bursters LMXBs V801 Ara (=4U 1636-536) and V926 Sco (=4U 
1735-444) through NIII A4640 fluorescent emission caused by ir- 
radiation. 

• The narrow NIII A4640 spots in the Doppler maps define 
K em = 277 ± 22 km s" 1 (V801 Ara) and K cm = 226 ± 22 
km s _1 (V926 Sco), and a new set of spectroscopic ephemerides 
which lend support to the assumption that photometric modula- 
tion is driven by the visibility of the irradiated donor stars. On the 
other hand, the phasing of the radial velocity curves of Hell A4684 
and the Bowen blend suggest that these lines are mainly associated 
with the disc bulge. Our spectroscopic ephemerides, combined with 
burst oscillations data for V801 Ara, enables us to refine the neu- 
tron star's projected velocity to K\ — 90 — 113 km s _1 for this 
particular system. 

• Following Muno z^Darias et all 120051) . we have computed the 
K — corrections for the two LMXBs and obtain K% = 360 ± 74 
km s _1 , q = 0.21 - 0.34, f(M) = 0.76 ± 0.47 M (V801 Ara) 
w&K 2 = 298±83kms" 1 ,g = 0.05-0.41, /(M) = 0.53±0.44 
Mq (V926 Sco). Both systems are seen at intermediate inclination 



angles in the range i ~ 30 — 60° , with V801 Ara having the higher 
inclination of the two. 

• Regarding V801 Ara, disc flaring angles a < 8° seem to be 
ruled out because of the high inclinations implied. Opening angles 
a ~ 12° support massive neutron stars and main-sequence donors 
which may descend from intermediate-mass X-ray binaries as pre- 
dicted by some evolutionary models (Schenker & Kins 2002). Al- 
ternatively, higher opening angles a ~ 16° are consistent with 
canonical neutron stars and main-sequence (or slightly evolved) 
donors which have evolved from standard LMXBs through a dy- 
namically stable mass tranfer phase. 

• The lack of an X-ray mass function in V926 Sco prevents to 
set tighter constraints to the K — correction, mass estimates and 
evolutionary history in this LMXB. 
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